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O ■ Abstract 

(^ • We write down the Yang-Mills partition function and the average Wilson loop in terms 

^^ . of local gauge-invariant variables being the six components of the metric tensor of dual 

space. The Wilson loop becomes the trace of the parallel transporter in curved space, 
else called the gravitational holonomy. We show that the external coordinates mapping 
the 3d curved space into a flat Qd space play the role of glueball fields, and there is a 
(D i natural mechanism for the mass gap generation. 
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1. The problem of formulating Yang-Mills (YM) theory in terms of gauge-invariant 
variables has been with us for the last 25 years. Several approaches have been suggested to 
trade the YM potential for some gauge-invariant variables. However, usually it is difficult 
to incorporate matter in such approaches as it couples not to gauge-invariant variables but 
to the gauge potential A^. 

The simplest but typical matter is the probe source, or the Wilson loop. It is given by a 
path-ordered exponent of A^, and it is difficult to present it in terms of any gauge-invariant 
variables. At some point it has been suggested that Wilson loops themselves are the only 
reasonable gauge-invariant variables, and the theory should be reformulated in terms of loop 
dynamics [jl], Q. Because Wilson loops are non-local objects there have been no decisive 
success on this way, in spite of enormous efforts taken in twenty years. 

In this Letter we suggest an approach which is opposite in spirit to loop dynamics. 
Namely, we reformulate the YM theory in terms of local gauge-invariant variables, such 
that Wilson loops are also expressible through these variables. We demonstrate it in a 
simple case of SU{2) theory in d = 3 Euclidean dimensions. Four dimensions and higher 
groups are, most probably, doable as well, and we comment on that at the end of the paper. 



2. The basis for our analysis is the reformulation of the SU{2) YM theory in d = 3 
in terms of quantum gravity. This has been performed both in the continuum p, ^, |^ by 
using the first-order formalism and on the lattice [^, ^ by making duality transformation on 
the lattice and then passing to the continuum limit. We start by briefly outlining the first 
derivation. 

The YM partition function can be written with the help of an additional Gaussian in- 
tegration over the 'dual field strength' G"; (a=l,2,3 are color and i=l,2,3 are Euclidean 
indices): 



Z = J DAt exp -^ J d'xFt^{A)Ft^{A) [Ft,{A) = d,A^^ - d.A^ + e'^^^A'lA 



DGt DA'^ exp fd^x 



.^G^G^ + ^^^'^F^^iA)Gl 



(1) 



Gaussian integration over A'^ can be easily performed, the saddle-point A" being defined 
by the equation 

^Vk^Q,^ab ^ ^acb^c-^^b^ ^ ^ijk D^b^A)Gi = 0. (2) 

This equation is generally satisfied if 

Df{A)Gl = r],Gt (3) 

where F is symmetric in subscripts. Eq. (|]) is known to be satisfied in general relativity, 
provided one makes the following identification: 

• G" = e^ where e" is a dreibein such that the metric tensor is the standard 

g^, = €'}€';, g^^ = e'^'e-^, e^e^^ = <5«^ det e^ = y^; (4) 

• r is the standard Christoffel symbol, 

^i,jk = - {dkQik + djgik — diQjk) , r*;, = 5f* Tijk] (5) 

• A'^ = — |e"'"^co'"'' where u is the spin connection, 

uf = -^^ = \e^'{d,el - dkc^) - ^e''{d,et - d,et) - ]f^ e''^ e\{d^e] - d,el). (6) 

It is easy to check that with the above definitions eq. (^ and hence eq. (H) are fulfilled. 
Eq. (|) gives an explicit solution for the saddle-point potential A" through the dual field 
strength G1 = e\. 

Pursuing further the analogy with Riemannian geometry one can define the covariant 
derivative in curved space, 

(v.)f = d^sf + rf„ (7) 

whose commutator is the Riemann tensor, 

[V.V,]f = i?'$,,. = d,T^,-d,Tl + TtTY,-T^^T^ (8) 

= e^''^Fr^{A)e'''e1, e? ^ G?, (9) 



related, as we see, to the YM field strength at the saddle point A. More directly, we get from 
eq. @ the YM field strength at the saddle point expressed in terms of gravity quantities, 

1 



F^M) 



ifu, ^ 



abc b cl 



(10) 



Therefore, the second term in eq. (|l|) is, at the saddle point, nothing but the Einstein-Hilbert 
action: 

e^^'F^^{A)Gl = ^R, R^R^^J\ (11) 

where we have used the identity 



abc a b 



y/g^kr. 



(12) 



The first term in eq. (||) is 

GtGt = e\el = gu- (13) 

3 3 

Integrating over A'} in eq. (|I|) gives a pre-exponential factor (detG'")~2 = g-i. The 
resulting integral in eq. (|l]) is over the 9 components of the dreibein G^ which can be rotated 
by the 3-parameter gauge transformation O""^. However, both the action and the measure 
depend on G^{x) only through the metric tensor Qij which is gauge-invariant. Writing 
^(9)ga _ (ii^)Q(^i> (i(^)g^jg-2 Que cau pass to integrating over the metric. One gets finally the 
YM partition function (|l|) identically rewritten in terms of 6 gauge- invariant variables 0, Q : 



Z = Dgij g 4 exp / d^x 



2 

-^gu+'-VgRig) 



(14) 



The second term, the Einstein-Hilbert action, is diffeomorphism-invariant and depends 
actually only on three variables. The first term is not, and we call it the 'ether' term ||^ as it 
corresponds to an isotropic and homogeneously distributed 'matter' which, however, spoils 
the general covariance. It distinguishes the YM theory from the pure gravity theory (which 
is a topological field theory [0]), and is responsible for the propagation of transverse gluons 
in the perturbative regime 0. We note that the integration measure also differs from the 
diffeomorphism-invariant measure Y[i<j dgijg~^. 

It is instructive to rewrite the partition function ([T^) further introducing 6 external 
coordinates w"{x), a = 1, ... 6, describing the embedding of a general 3-dim Riemannian 
manifold into 6-dim flat space. The metric induced by embedding is 



gij = diW^djW^ 



(15) 



Quantities F*^, y^ and R written in terms of w" can be found in ref.^. 

The Jacobian for the change of integration variables from the 6 components gij{x) to 
6 external coordinates w°'{x) is the 6x6 determinant made of second derivatives, Dg^j = 
Dw;"det(w"). It is zero in the degenerate case when the manifold can be embedded in 
less than 6 dimensions. It terms of the external coordinates the YM partition function has 
the following form (it is in this form the partition function has been derived in ref.0 from 
making a duality transformation on the lattice and then passing to the continuum limit): 



Z= / Diu'^ix) det{didjw'') g( 



w) 



exp 



d^ 



X 



2 

■^d,w''d,w''+'-^/^)R{w) 



(16) 



3. Let us now consider the average of the Wilson loop in representation J = |, 1, |, 
along a closed contour a;*(r), defined as 



Wj 



2J + 



1 f Ht^ 

Tr P exp i it/r — A"" T" 

+ 1 ^ J dr ' 



(17) 



where T"^ are the SU{2) generators in representation J. It is not at all clear that one can 
rewrite it in terms of local gauge- invariant variables Qij or w'^. A hint that it may be possible 
comes from the dual formulation on the lattice where the Wilson loop appears to be a product 
of certain 9j symbols along the contour 0, each of which is gauge- invariant. We shall show 
that in the gravity formulation the average Wj becomes the average of a parallel transporter 
in curved space, else called the holonomy, taken in the appropriate representation (spinor, 
vector, etc.) 

To that end it is convenient to write down the Wilson loop in a form where the path 
ordering is traded for an additional functional integration over a unit 3-vector n 'living' on 
a surface spanned on the contour [H, 



W, 



Dn(a, t) exp 



iJ 



d^%i^>'') + Y 



d'S'^ e^'^n^dm'dm' 



Dn(a, t) 



exp 



- / d^S'^ 



{-F^^^n'' + €'"'' n^{Dinf{D,n) 



'M 



(19) 



This representation for the Wilson loop is convenient because A" comes linearly in the 
exponent. The second line is a 'non-Abelian Stokes theorem'; it is written down in an 
explicitly gauge-invariant form but actually the terms quadratic in A"} cancel out. 

One defines the average Wilson loop by inserting eq. (|l^) or, equivalently, eq. (|T8D into 
the partition function (|I]). Since the integral over Af is Gaussian and since the source term 
Afn°' is linear in A"}, it is sufficient to replace Af in the Wilson loop by the saddle-point 
value Af. One can disregard the quadratic point-like self- interacting term emerging from 
this replacement, as being independent of the dynamics. 

Meanwhile, Af is given by the spin connection, see eq. (^. The path-ordered exponent 
of the spin connection determines the parallel transporter of spinors in curved space. Par- 
allel transporters of higher spins can be always built from spinors. We, thus, come to the 
conclusion that in the quantum-gravity formulation of the YM theory one can calculate the 
Wilson loop as the trace of a parallel transporter along a closed contour in curved space. 

4. Spin connection u;°* is defined by the dreibein ef which is not uniquely determined by 
the metric tensor. In principle, for a given metric one can use any dreibein to calculate the 
Wilson loop, but this is not very satisfactory. Since Wilson loops are gauge-invariant and so 
is the metric tensor one would expect that parallel transporters are expressible through the 
metric tensor (and its derivatives) only. 

Let us first consider Wilson loops in integer- J representations. In this case the problem 
of expressing them through the metric tensor is immediatelly solved by the observation |^ 
that the holonomy defined via the spin connection is equal to that defined via Christoffel 
symbols. For example, for the J = 1 (vector) representation one has |^: 



w 



J=l 



exp 



, dx^ 

dT-—UJi 

dr 



/dx^ 
dT — Ti{g) 



(20) 



For half-integer representations the problem is more difficult as there is no way to define 
the parallel transporter of spinors other than via the spin connection. Nevertheless, this 
problem has been solved in ref.||^ with the help of the 'gravitational non-Abelian Stokes 
theorem' similar to eq. ([T9|). It is only for parallel transporters of spinors along a closed loop 
and only when written down in a surface form that the spinor holonomy can be expressed 
through the metric tensor and its derivatives. Actually, one can write a unified formula for 
the holonomy for any spin J; it is given by a functional integral over a covariantly unit vector 
m 'living' on a surface spanned on the contour [Q: 

Wj = J Dm y/g b{m'vn? gij - 1) 

X exp t^ JdSk ^ [{R6^ - 2R';) m^ + e''" e„r ^^(V. m)«(V, my] (21) 



where R^ is the Ricci tensor, R = R\is the scalar curvature and Vj is the covariant derivative 
in curved space (|^. Despite its surface form, eq. (pi]) is independed of the way one draws 
the surface. 

Eq. ( pT] ) gives a general solution to the problem of expressing the YM Wilson loop in terms 



of the metric tensor; it should be inserted into the partition function ([1^) for calculating the 
average Wilson loop over the ensemble of YM configurations. If external coordinates w'^{x) 
are used as dynamical variables instead of gij {x) one has to express R\. and F* ^ through w" 
first, and plug them in eq. (|T^. In both cases the average Wilson loop is expressed through 
local gauge- invariant variables, be it Qij or w'^. 

5. So far we have presented exact results, and now we would like to speculate on how the 
expected mass gap for glueballs and the area behavior of the Wilson loop might be achieved. 
We also give an example where eq. ( PT| ) is considerably simlified. 

A characteristic feature of the partition function ([T6|) is the purely imaginary Newton 
constant. It damps field configurations leading to fast oscillations. A trivial way to get rid 
of those oscillations is to put i? = 0, which leads to the perturbative regime p. A less trivial 
example of how to avoid fast oscillations is to have the Einstein-Hilbert action quantized 
so that the oscillating factor becomes exp(27rm) = 1. This is achieved e.g. in de Sitter 5*^ 
spaces of constant (quantized) curvature R = 547r^/n^. 

Indeed, the de Sitter space with curvature R can be generally parametrized by 4 external 
coordinates w"^{x)^ A = 1, . . . , 4, restricted to the S^ sphere of radius JQ/R: 

J: w\x)w\x) = -|. (22) 

A=l ^ 

The square root of the metric determinant becomes a determinant itself being equal to 

^g = iMe^^' 6^^^^ d,w^ d,w^ d.w'^ w^, (23) 



so that the Einstein-Hilbert term in eq. (p!6| ) is 

'- I d'x^R = I GTr'J^N^i^d, (24) 

Z J V -K 



where iV„ind is the properly normahzed winding number of the mapping 5''^ ^—>- S^ reahzed 
by w^{x). Putting it to be an integer n and demanding that the l.h.s. of eq. (0) is 27rm 
we arrive to the quantization condition written above for the scalar curvature. Such values 
of R are preferred because there are no damping oscillations then. 

The 'ether' term is the kinetic energy term for four fields w'^^x) subject to the constraint 
(p2D, therefore we are left with a low-energy 0(4) sigma model. The large- iV sigma models 
are known to develop a mass gap; technically it is achieved by introducing a Lagrange 
multiplier for the constraint and then finding a saddle-point value of this multiplier, which 
is justified at large A^. It is not clear if one can consider 4 as being much much larger than 1 
but this scenario seems probable. That would give a mass gap to glueballs which we identify 
with the w"^ field. 

Whether this dynamical regime leads to an area law for the holonomy in half-integer 
representations, is less clear. On one hand it is intuitively a rather natural outcome: the 
constant curvature is a gauge-invariant analog of a constant field strength usually leading to 
the area law. On the other hand we have not been able to prove it directly so far, and the 
question remains open. 

We would like to remark that a const ant- curvature background is a gauge-invariant analog 
of the long- sought master field leading to a nonzero gluon condensate. 

Another interesting regime is a cylindric space S"^ x TZ also of constant curvature R but 
with one eigenvalue of the Ricci tensor being zero. Such spaces are generally parametrized 
by four coordinates w with only three of them subject to the constraint 



a=l 



w^Y = — , w'^ arbitrary. 

R 



(25) 



In this case we have managed to find a general expression for an arbitrary Wilson loop. It 
is given by the SU{2) character whose argument is the winding number of w"'{x) M: 



Wj 



2 J 



J 

E 

m=—J 



"Y E exp 



im^. 



$ 



2_Y 1 
r) 2 



J dSk eabc e'^' d.w'' d, 



w^w\ 



(26) 



where dSk is the dual element of the surface spanned on the contour. We note that eq. (|26[ 
gives an example of the Wilson loop being expressed through local gauge-invariant variables 
though it is less general than eq. (pTf ). 

For cylindric spaces the Einstein-Hilbert action is a full derivative. 



d^Xy^R 



-R2 



d xdh 



^4^iifc^abc^.^a^^.^b^. 



(27) 



and can be reasonably set to be zero, therefore there are no oscillations in the partition 
function, too. The 0(3) sigma model we are left with in this case may also develop a mass 
gap for glueballs. However, this regime most probably does not lead to the area law. The 
reason is that for large areas inside the Wilson loop the fiux $ in eq. (^) becomes the 
winding number of the mapping S"^ i— >■ S*^ but with the coefficient 47r, therefore non-trivial 
mappings w°'{x) are not 'felt' even in half-integer representations. A more careful analysis in 
the spirit of the CP^~^ model shows that the static potential has a logarithmic, not linear, 
behavior. 



6. In conclusion, we have reformulated both the YM partition function and the average 
Wilson loop in terms of local gauge-invariant variables, being actually the metric of the 
dual space. It looks as a promising starting point for studying two features expected in YM 
theory: the mass gap for glueballs and the area law. Generalization to higher gauge groups 



and to d = 4 is possible e.g. along the lines suggested in refs.M, 11 1. In rf = 4 a formula 



similar to eq. (|2l|) exists, expressing the holonomy through the metric tensor 
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